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Abstract 

A nonperturbative calculation of elastic electroproduction of the J/tp meson is presented and 
compared to the experimental data. Our model describes well the observed dependences of the 
cross sections on the photon virtuality Q 2 and on the energy, and the measured ratio R = ctlI^t 
of longitudinal to transverse cross sections. 
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I. INTRODUCTION 



In a previous paper |l[ we have investigated photoproduction of heavy mesons off protons, 
finding good agreement between experiment and our calculations based on a nonperturbative 
approach to high energy scattering 0,0], where no free parameters had to be introduced. 
In this letter we present the results for electroproduction of J/ip mesons, and compare our 
results with the published HERA data. 

In the perturbative approach to electroproduction of heavy mesons the coupling of the 
exchanged gluons to the heavy vector mesons is treated perturbatively, while for the coupling 
of the gluons to the proton an external nonperturbative quantity, the gluon density in the 
proton, has to be introduced. In our approach to photo- and electroproduction of heavy 
mesons, the small system, the 7* — J/tj) transition overlap, and the large system, the proton, 
are treated on the same footing and the treatment of high energy scattering is based on 
functional integrals 0, 0], which are evaluated in a nonperturbative approach to QCD, the 
stochastic vacuum model[a,E3|- It has been shown 0] that this nonperturbative approach also 
implies factorisation in the sense that genuine nonperturbative effects in the small system 
can be absorbed into the large system, that is the proton in the present case. Although 
our nonperturbative approach is more model dependent than the perturbative calculation, 
it offers the advantage that purely hadronic reactions are described with the same set of 
parameters as photo and electroproduction processes. Therefore in our calculations of heavy 
vector meson production no new parameters are introduced and the Regge model can be 
applied directly. Furthermore the influence of confinement effects in the small system can be 
studied. The WKB framework underlying our model also allows to calculate the dependence 
of the production processes on (moderate) momentum transfer, which is not directly possible 
in the perturbative approach. 



II. BASIC FORMULA AND GENERAL RESULTS 

For convenience we present here some basic formulae developed in our previous work on 
photoproduction where motivation and details can be found. 
The electroproduction amplitude of vector mesons is written 

Ty* p -+v P ,x{t) = J d 2 R l dzii) V x{zi,Ri)*i) 1 *x{z 1) R l ,Q 2 )J(q,z 1 ,R l ) , (2.1) 



with 



J(g,zi,4) = / d 2 R 2 d 2 be^%(R 2 )\ 2 S(b,z 1 ,R 1 ,z 2 = 1/2,4) • (2.2) 



Here S(b, zi,R\,l/2, R 2 ) is the scattering amplitude of two dipoles with separation vectors 
Ri, R 2 , colliding with impact parameter vector b; q is the momentum transfer of the reaction 

t = -q 2 - m J(Q 2 + M 2 V ) js 2 + 0(s- 3 ) « -q 2 . (2.3) 

The differential cross section is given by 

do 1 

\T\ 2 . (2.4) 



d\t\ 16its< 

The wave functions of the photon and vector meson have been discussed extensively in 
1]), where it has been shown that the results obtained with different forms of the meson 
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wave function are very similar. We use here only the Bauer- Stech-Wirbel (BSW) |8j type of 
wave function, which is of the general form 

ipvx(z\, Ri) — f( z i) exp[— uj 2 R\/2] x helicity dependent factors , (2.5) 

where f{z\) contains the normalisation constant and the 24-dependence of the meson wave 
function, with the form 



N 



f(z) = —^/z(l-z) exp -^L( z --y . (2.6) 



The size of the vector meson is determined by u, which is fixed by the leptonic decay 
width of the meson. 

After summation over helicity indices, the overlaps of the photon and vector meson wave 
functions are given by 



2^ 



/Vy,±i(^i,#i,<2 2 ) = Vv,±i(zi,^i)*Vy,±i(2i,-Ri,(9 

271 

x (e ^R^zl + (1 - ZxfjK^e R x ) + m 2 f K (e R x ] 
Pr,v,o(zi, Ri, Q 2 ) = ifo,o{zi, R 1 )*ip r>0 (z 1 , R u Q 2 ) = (2.7) 



-16 u f(z 1 )eM-LU 2 R 2 1 /2]z 2 (l-z 1 ) 2 QK (e R x ) 

271 



Here 



I Zl (l - Zl )Q 2 + m 2 , 

Ko, K\ are the modified Bessel functions, A = ±1 and denote transverse and longitudinal 
polarisations of the vector meson and the photon, is the quark mass and e/ is the quark 
charge in units of the elementary charge for each flavour /. All parameters are fixed from 
other processes (see fH) and all observables can be calculated from (j2.1j) . 

The energy dependence in our model is motivated by the two-pomeron model of Don- 
nachie and Landshoff 0. For R\ < r c w 0.22 fm the coupling through the hard pomeron 
induces the energy dependence (s/s ) 0A2 , while the coupling of large dipoles follows the 
soft pomeron energy dependence (s/sq) 0,0808 . The reference energy is So = (20 GeV) 2 . We 
therefore split the integration over Ri appearing in (J2.1)) into two parts, as fully described 
in our study of J /if) photoproduction 

Before proceeding to the comparison with experimental data we wish to present some 
general approximate features of our model which describe the overall situation quite well 
and provide a background against which finer details can be studied. 

If Ri is small compared to the extension of the proton and \t\ = q 2 < 1 GeV 2 , we obtain 
after integration over the azimuthal angle of R\, the simple expression for the amplitudes 

T^ Vp {t) ~ (-2*5)2^^-—- J dR x dz x Rle-^p r y, x ( Zl , R h Q 2 ) . (2.8) 

For the case of J / if) production the numerical values for the constants in ()2.8|) are 

C = 2.28, a = 0.847 GeV" 2 , b = 3.85 GeV" 2 , rj = 0.059 . (2.9) 
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This result means that the integrated expression in (|2.2|) is proportional to 
R\ exp[— ?7-Ri|t|] and that the main t dependence in the amplitude is factorized out, be- 
ing independent of Q 2 and of W. 

Since in the overlap functions with heavy vector mesons i]R 2 is very small , the exponential 
factor e~ vR ^ can be extracted from the integral in ([2.8)1 . and represented by an external 
factor e _r?i?m '*', where R m is an appropriate mean value of R\. The numerical calculation 
shows that a very good representation for the Q 2 dependence of this external factor is given 
bjR 2 m ^12/(Q 2 + M 2 ) . 

We can then write for the scattering amplitude 

-12 V \t\/(Q 2 + M 2 ) 

x ( ThX{Q2 Ki) eh +TsX(Q2 KiT) (2 - 10) 

with 

T£(Q 2 ) = 2tt £dR 1 dz 1 (^pj Rlp JtVA (Q 2 ,z 1 ,R 1 ) , 

POO 

T S \Q 2 ) = 2tt / dRxdzxRlp^v^.zuRx) . (2.11) 

J r c 

The parameters for the energy dependence are s = (20 GeV) 2 , e h = 0.42, e s = 0.0808 . 

Eq. ()2.10|) allows us to express all observables through the functions T^{Q 2 ) and T^{Q 2 ) 
(|2.1ip . We display these amplitudes in Fig. ^ They can be parametrized for the range 
< Q 2 < 100 GeV 2 in the forms 

T±\Q 2 ) 

t£ 1 (Q 2 ) 



T^ Vp (s,t;Q 2 ) w (-2is)— -— exp 



0.0267/(1 
0.00488/(1 + 



Q 2 



\2.71 



15.03' 

Q 2 

20.46' 



vl-72 



(2.12) 



T°M) 

T° h (Q 2 ) 



-0.00635Q/(1 + 
-0.00170g/(l + 



Q 2 



\3.23 



18.00' 

Q 2 , 
20.61' 



1.9 



where Q 2 is measured in GeV 2 . 

The overall \t\ dependence determined by the form factor 



F(\t\) 



-a\t\ 



l + b\t\ 



exp[-12^|/(Q 2 + M 2 )] 



(2.13) 



yields a characteristic curvature in the log plot. We can introduce an effective slope defined 

as 



B 



da 



d\t\ 



o 



/ d\t\ 



da 

d\t\ 



(2.14) 



which is the inverse of the integral / -F(|t|) 2 (i|t| and varies monotonously from B 
7.038 GeV -2 to B = 6.839 GeV" 2 in the interval < Q 2 < oo. 
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T - transverse 
L - longitudinal 



Q 2 (GeV 2 ) 
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FIG. 1: Hard and soft parts of the reduced amplitudes for transverse and longitudinal polarisations, 
according to Eqs. Q2.11JI . These functions can be represented by the parametrisations (|2.12j) . (|2.13j) . 



In any approach where the wave function of the heavy meson is taken into account 1(1 11 



nonperturbative quantities like the "meson radius" 1/u enter even at the "small" side of 
the interaction which determines the "hard" scale. It is therefore interesting to study the 
strictly non-relativistic limit of the meson wave function, where the ratio uj/My goes to zero 
and the meson radius drops out in the final expressions for the amplitudes. In this limit, 
which also implies rrif — > My/2, all integrals can be performed analytically if we neglect the 
energy dependence by putting e s = €h = 0. We obtain the asymptotic results 



T^ Vp (t,Q 2 )/(-2 lS ) 



Ce- a ^ V3 



16 



1 + 6|t| y/a {Q 2 + M£) 2 



M 3/2 r l/2 
1V1 V 1 e+e- J 



"f*p—>Vp 



(t,Q 2 



Q 



y Jt, Q 2 



M -rp-+v P \ 



In Fig. |21 we show in solid line the full result (with = e s 
asymptotic form of ()2.15|) for forward differential cross section 



(2.15) 

0) and in dashed line the 



da/d\t\ 



(t = 0). 



III. NUMERICAL RESULTS AND COMPARISON WITH EXPERIMENTS 

Our numerical calculations presented in the figures below with experimental data 
are exact evaluations of the model, based on Eqs. (j2.1l2.2|) . However, we stress that 
Eqs. (12. 1012.1112. 1212. 13|) represent very accurately these results. 
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FIG. 2: Influence of the nonperturbative parameters of the J/ip meson on the Q 2 dependence 
of the forward differential cross section. In solid line our calculations without energy dependence 
(e^ = e s = 0) and in dashed line the asymptotic form of Eq. (|2.15|) . 

In Fig. |3]we show the integrated J/ip production cross section o = ctt+ctl at W = 90 GeV 
as a function of Q 2 , together with the HERA data from the HI and Zeus collaborations 
[I2I fl3| . The solid line corresponds to the exact equations (j2.H2.2J) . while for the dashed 
line there is a multiplying factor 

4tt/0.57 
lllog(Q 2 + 7.42) 

introduced [3] to take into account the running of the strong coupling in the model which 
otherwise is purely nonperturbative. 

The energy dependence adopted in our model is fully compatible with the existing pub- 
lished data jia . Eil . EH [Tril l as can be seen in FigEJ We repeat in this figure our results for 
photoproduction (empty circles and squares) which have been shown before [1] to agree well 
with the data. 

The energy dependences of the cross section have general forms 

*tAQ 2 ) = K£ d (Q 2 ) W 2 ^ + A T J t {Q 2 ) W 2ea ) 2 , 

where A^ rd (Q 2 ) and A^ t (Q 2 ) can be easily calculated from (ETTUJ) . (gUD and (|2~T31) . 

Often experimental cross sections are fitted through a single power 
which must be used in a limited energy range. For small Q 2 the soft pomeron is important 
and 5 is small, but for larger Q 2 the hard pomeron becomes dominant and the exponent 
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y* p -» y p 




W=90 GeV 



squares - 1-11(1999) 
circles -ZEUS(1999) 

solid - theoretical calculation 



dashed - with a(s) factor 
(4*3.1 4/0.57)/(11*LOG(x+7.42)) 
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FIG. 3: Q 2 dependence of the integrated elastic cross section a = ot + <Jl a t the energy W = 
90 GeV, compared to data from Zeus 0] (circles) and HI 0] (squares). The full line represents 
our results as explained in the text. In the dashed line an extra factor is introduced to account for 
the possible effect of a Q 2 dependence of the running strong coupling constant [14 1 . 



5 = 4e c g — ► 4eh = 0.168 as Q 2 — > oo. A convenient interpolation formula in the range 
20 < W < 200 GeV is ^ 

e cS (Q 2 ) = e h - {3 28 + Q2/M 2y.22 ■ 

Our results for the transverse and longitudinal cross sections at W = 90 GeV can be 
parametrised in the forms 



63 50 Q 2 / M 2 



The ratio R = <Jl/<Jt of longitudinal to transverse cross section is displayed in Fig. El The 
solid line shows our result for W = 90 GeV, which agrees well with the data [l2l . at the 
same energy. The dashed line is the asymptotic result -R a sym P = Q 2 /^ 2 ^, from (|2.15|) . 

The differential cross section da/d\t\ and its energy dependence obtained from our model 
have been shown ^ to describe well the data for J /if) photoproduction in a wide energy 
range. There are no published data of da/d\t\ for electroproduction to be compared to our 
calculations, which predict that the shape of the angular distribution is almost independent 
ofQ 2 . 
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W (GeV)° 



FIG. 4: Energy dependence of the integrated cross section for several values of Q 2 compared to 
data from Zeus in full ^3] an d empty circles and from HI in full and empty squares. 



IV. CONCLUSIONS 



Our nonperturbative model describes the data well. We found that the main features 
are expressed in the simplified equations (|2. 11)12.1112. 1212.1"!?)) which represent very well the 
full results of our model. This is due to the small range of the overlap functions compared 
to the proton size and to the range of the correlation functions of the QCD vacuum , 
which determine the interactions in our model. As a consequence of these properties of 
the amplitudes, in our nonperturbative treatment of the electroproduction of heavy vector 
mesons (here with particular application to J/if> production) some quite general features 
emerge, which are mainly a consequence of the general approach to high energy scattering 
0, IE HI and not of the specific stochastic vacuum model 0,0], which yields the numerical 
values of the parameters in (|2.9|) . 

The main shape of the Q 2 dependence is determined by the overall overlap strength (|2.11|) . 
Our calculations, which reproduce well the data for J /if} production for Q 2 < 60 GeV 2 devi- 
ate considerably from the asymptotic result (|2.15|) and shows that genuine nonperturbative 
parameters, like the size of the vector meson (1/u), have a considerable influence even on 
the "hard part" of the interaction. 

Fig. H| shows that the PF-dependence based on the two-pomeron model || is appropriate 
to describe the existing data. 

In the ratio R = o^j 'a? the characteristic numbers coming from the stochastic vacuum 
model nearly cancel out, as can be seen in (|2.8|) . and the ratio depends solely on the overlap 
strengths defined in (|2.11|) . Due to the energy dependence induced by the overlap function, 
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FIG. 5: Ratio between longitudinal and transverse cross sections at W = 90 GeV as a function of 
Q 2 compared to the data from Zeus (circle) and HI 0] (squares) collaborations. The solid 
line is our theoretical result, for which we give a parametrization. The dashed line represents the 
asymptotic calculation (size zero for the vector meson). 

see (J2.1UI2.11|) . the ratio will also exhibit a characteristic energy dependence. At W = 
90 GeV the ratio agrees reasonably well with the data shown in Fig. It is important to 
check this result against future data at other energies and for a wider range of Q 2 . 

The predicted t dependence is nearly universal, namely independent of Q 2 and W, with 
only small corrections due to the term e _r?Rl '*' in f)2.8j) and slight dependence on W due to 
the different energy factors controlling the couplings of small and large dipoles. 
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